Stripes, Clusters, and Nonequilibrium Ordering for Bidisperse Colloids with 
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We show that two-dimensional bidisperse assemblies of colloids with strictly repulsive interactions 
exhibit stripe, cluster, and partially crystallized states when driven over a quenched random sub- 
strate. The nonequilibrium states on a substrate are significantly more ordered than equilibrium 
states both with and without substrates. A minimum substrate strength is necessary to induce the 
nonequilibrium pattern formation. Our results suggest that a combination of driving and quenched 
disorder offers a new approach to controlling pattern formation in colloid mixtures. 
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Understanding how assemblies of particles organize is 
an outstanding problem in both equilibrium and nonequi- 
librium systems. The ability to control pattern formation 
would have a profound impact on applications that use 
large-scale self-assembly to create specific pattern mor- 
phologies. In equilibrium, it is known that particles with 
competing repulsive and attractive interactions can or- 
ganize into clusters, stripes, and spongelike textures 
Similar patterns can also arise for particles with strictly 
repulsive potentials that have a two step form Re- 
cently, it was shown in both simulations and experiments 
that a two-dimensional (2D) bidisperse assembly of col- 
loids interacting via repulsive magnetic dipoles of two dif- 
ferent strengths can cluster into spongelike patterns [|[. 
This suggests that other types of ordering are possible 
for bidisperse repulsively interacting particles. 

Here, we demonstrate that a variety of distinct par- 
tially ordered states can occur for a bidisperse system of 
strictly repulsively interacting colloids driven over a ran- 
dom substrate that is sufficiently strong. There has been 
extensive work on the ordering of monodisperse charged 
particles moving over random quenched disorder for sys- 
tems such as vortices in type-II superconductors 0, H, Q , 
moving Wigner crystals 0] and colloids In the ab- 
sence of quenched disorder, a monodisperse system forms 
a triangular lattice, while strong quenched disorder dis- 
torts the lattice and generates numerous topological de- 
fects. These defects can be partially annihilated and 
the order partially regained by driving the system into 
a nonequilibrium moving smectic state [!, @, 0, Hi- 
lt is not known what type of nonequilibrium configura- 
tions would occur for two species of particles driven with 
identical force over quenched disorder. In contrast to the 
monodisperse case, for a bidisperse system even the equi- 
librium states in the absence of a random substrate are 
intrinsically disordered Q , so one might expect the addi- 
tion of quenched disorder to further disorder the system. 
Instead, we find that for certain regimes, nonequilibrium 
bidisperse colloid assemblies moving over quenched disor- 
der have more topological order than the corresponding 
equilibrium states, even those without quenched disorder. 
The quenched disorder must be sufficiently strong for the 



topological ordering to occur. We argue that for strong 
quenched disorder, significant plastic deformations of the 
driven colloid configuration occur due to the fact that the 
different species move at different average velocities close 
to the depinning threshold. We show that the system 
can organize into stripes with triangular ordering within 
each stripe. When the difference between the two species 
is small, a moving smectic state forms, whereas when the 
disparity is very large moving cluster states form. Our 
results are robust for a wide range of colloid densities, 
quenched disorder strengths, and system sizes. 

We employ Brownian dynamics to simulate a 2D sys- 
tem of size L x L containing N colloids and N p poten- 
tial traps with periodic boundary conditions in the x 
and y directions. The overall density of the system is 
n = N/L 2 . The colloids interact via a screened Coulomb 
potential V(Rij) — (Eq / Rij) exp(— nRij) , where R,*(j) 
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is the position of colloid Rij 

Z* 2 /(Aneenan), Z* is the unit of charge, e is the solvent 
dielectric constant, and 1/k is the screening length. The 
unit of distance in the simulation is do, and unless other- 
wise noted, L = 48ao. In this work we fix k = 4ao, which 
is reasonable for experiments on colloids in nonpolar flu- 
ids @ . Forces are measured in units of Fq = En. /an. The 
dynamics of colloid i are given by the equation of motion 
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where rj is the damping term, q^j\ is the dimensionless 
charge of colloid Rik = |R» — Rfc|, and Rk is the 
position of trap k. In this model, hydrodynamic effects 
are neglected; such effects can be strongly screened in 
a system confined within 2D walls. To introduce bidis- 
persity to the system, half of the colloids have charge 
Qi = a a and the other half have qi — qs, where we fix 
qs = 1. The quenched random substrate is modeled by 
randomly distributed parabolic traps of density n p and 
radius r p = 0.1a with V p (R ik ) = ~(F p /2r p )(R ik ~ r p ) 2 
for Rik < r p and zero interaction for Rik > r p . Here F p 
is the maximum pinning force. This model for quenched 
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FIG. 2: (a) Average velocity (14) vs applied drive Fd for the 
system in Fig. 1. Light upper line: sample with F v — 0; dark 
lower line: sample with F v — 1.0; dashed line: d{V x )/dFr> for 
the F p — 1.0 case, (b) Pe vs Fd for the same system. Upper 
curve: F v = 1.0; lower curve: F v = 0. 



FIG. 1: A bidisperse system of colloids with qa/^b = 3 
and n = 0.385/(Jq. (a,c) Real space colloid configurations. 
Small circles: species A; large circles: species B. (b,d) The 
corresponding structure factor S(k). (a,b) Equilibrium state 
with no pinning or driving. (c,d) Moving stripe state for F v — 
1.0, n p = 4.0/oq, and F D /F C = 3.0. 



random disorder has given results comparable with other 
models in monodisperse particle systems. The externally 
applied driving force is identical for all particles and is 
given by F^j = Fdx. A uniform drive of this type could 
be created electrophoretically Thermal effects are 

modeled by random Langevin kicks with (F^f) — and 
(F?(t)F?{t r )) = 2r]k B T5 ij S(t - t'). The initial colloid 
configurations are obtained using simulated annealing. 
We then set F T = and gradually increase Fd in incre- 
ments of 5Fd — 0.001 every r = 5000 simulation time 
steps. After each drive increment, once the system has 
reached a stationary state we measure the average veloc- 
ity (V x ) = ((1/N) J2iLi v « ' where Vj is the velocity 
of colloid i. The depinning threshold Fq corresponds to 
the value of Fd at which (V x ) = 0.04Fo- As an example, 
for a = 0.6/iTO, e = 2, and Z* = 300e, F = 27.8 pN. 

We first consider a system where the ratio between the 
charges of the two colloid species is qA/lB = 3. In the ab- 
sence of quenched disorder, the colloids form a disordered 
mixed assembly after annealing, illustrated in Fig. HJa). 
To characterize the configuration, we determine the 

2 

structure factor, 5(k) = (1/N) Yli=i ex P( — *k ' Hi) ■ 
In Fig. [lib), S(k) for the equilibrium case with N = 864 
colloids at density n = 0.385/oq shows a ring struc- 
ture which is characteristic of disordered systems. We 
also analyze the fraction of six-fold coordinated col- 
loids Pq obtained from a Voronoi construction, Pq — 



(W)£-=i%i-6), where Zi is the coordination num- 
ber of colloid i. In a triangular lattice all the colloids 
have Zi — 6, giving Pq = 1.0. For the configuration 
shown in Fig. [TJa), Pq = 0.41, indicating that a large 
fraction of the colloids have Zi ^ 6. If we anneal in the 
presence of quenched disorder, we find a similar disor- 
dered state as measured by 5(k) and Pq. In general, for 
stronger quenched disorder the system becomes more dis- 
ordered, producing increased smearing in the ring struc- 
ture of 5(k) and reducing P 6 to P 6 = 0.3. 

In Fig. [DJa) we plot the dynamic response of the sys- 
tem given by the average colloid velocity (V x ) versus ex- 
ternal drive Fd for F p — (light line), and for a system 
with traps of strength F p — 1.0 and density n p = 4.0/ag 
(heavy line). For F p = 0, the velocity response is strictly 
linear. In the presence of traps, there is a clear depinning 
threshold Fq = 0.7 for colloid motion followed by a non- 
linear regime which crosses over to an ohmic regime at 
higher F D . For F p = 1, Fig. C^a) shows that d(V x )/dF D 
has a peak in the nonlinear regime at Fd = 0.75 and 
flattens in the ohmic region for Fd > 2.0. These features 
in the velocity-force curves are the same as those previ- 
ously observed for monodisperse systems with quenched 
disorder 0, 0, [|[ , where the onset of the ohmic response 
was correlated with a dynamical reordering of the system 
into a moving smectic state. 

In Fig. [DJb) we plot the corresponding Pq versus Fd 
curves for the systems with and without quenched disor- 
der. For F p = 0, Pq = 0.41 for all F D . For F p = 1.0, we 
still find Pq — 0.41 in the pinned region at Fd/Fc < 1. 
Above the depinning transition, Pq increases over the 
same range of Fd where there is a nonlinear response in 
the velocity force curves, as seen in Fig. [U^a). Here, the 
system moves plastically with portions of the colloids re- 
maining temporarily trapped while other colloids move 
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FIG. 3: (a) The different regimes for the system in Fig. 1 at 
n = 0.385/aj) and n p = 4.0/clq as a function of Fd and F p . I: 
pinned regime; II: moving disordered regime; III: plastic flow 
regime; IV: moving stripe regime, as in Fig. 1(c). Circles: the 
depinning threshold Fc ■ Squares: stripe formation threshold, 
(b) The different regimes for Fd vs density n at fixed F p — 1.0 
and n p — 4. 0/0,9. Symbols are the same as in (a), (c) Pa 
versus F p for n = 0.385/ao and n p = 4.0/oq at Fd/Fc = 3.0. 
Inset: Pe versus pinning density n p for fixed F p = 1.0 and 
n = 0.385/ag at F D /F C = 3.0. (d) Pe versus n for fixed 
F p — 1.0 and n p = 4.0/ao at Fd/Fc = 3.0. Inset: Pe versus 
system size L for fixed F p = 1.0, n = 0.385/ag, n p = 4.0/ag, 
and Fd/Fc = 3.0. 

around the trapped colloids. For Fd > 2.0, Pe saturates 
at Pq ps 0.7. This result indicates that the moving system 
with quenched disorder is more ordered than the equilib- 
rium system without quenched disorder. To illustrate the 
nature of this order, in Fig.QJc) we show the colloid con- 
figuration for the system with F p = 1.0 at Fd/Fc = 3.0, 
and we plot the corresponding S(k) in Fig. HJd). Here 
a stripe ordering occurs where the species have partially 
segregated. Complete segregation is prevented since this 
would produce a strong charge inhomogeneity in the sys- 
tem. Local triangular order appears within each stripe. 
The structure factor has a smecticlike form, indicative 
of the presence of stripes, as well as second order peaks 
at smaller scales resulting from the in-stripe triangular 
ordering. 

We next show that sufficiently strong quenched disor- 
der must be present for the stripe patterns to occur. In 
Fig. [3fa) we outline the four distinct regimes that we find 
as a function of driving force Fjj and pinning force F p for 
a system with n — 0.385/ap and n p — 4.0/oq. At all F p , 
the colloids are pinned in region I for Fd < Fc- As Fd 
is increased above Fc, we find that for F p < 0.4, the 
system never orders into a stripe state. Instead, the dis- 
ordered mixture depins with little or no plastic deforma- 
tion, and the moving state, marked region II in Fig.[3fa), 
retains the same disordered configuration as the pinned 
state. We have tested different values of SFd and find 
no change in the location of region II. For F p > 0.4, a 
plastic flow state (region III) appears above depinning. 



Here the colloids are repeatedly trapped and escape from 
the traps but stripes do not completely form. Finally, 
for high enough values of F p and Fd we find region IV, 
the moving stripe state. In Fig. |3{c) we plot Pe as a 
function of F p for the same system in Fig. [3ta) at fixed 
Fd/Fc = 3.0. For F p > 0.4, P e saturates near 0.7 and 
the moving stripe state forms. In contrast, for F p < 0.4 
the moving state has Pq = 0.41, the same value measured 
in the pinned state. This indicates that there is a critical 
disorder strength that is necessary for the formation of 
the stripe state to occur. 

In Fig. [3£b) we plot the locations of the four regimes 
as a function of Fd versus colloid density n for fixed 
F p = 1.0 and n p = 4.0/oq. The depinning threshold Fc 
marking the end of the pinned region I drops to lower 
values of Fd as n increases since the colloid- colloid in- 
teractions become stronger relative to F p , making the 
traps less effective at pinning the colloids. A reordering 
into the moving stripe state (region IV) occurs only for 
n > 0.2/oq; for smaller values of n, we find the moving 
disordered region II at high drives. This is confirmed by 
Fig. [3jd) , which shows Pq versus n for the same system 
in Fig. Mb) at fixed F D /F C = 3.0. For n < 0.2/a% in 
the moving region, P§ ~ 0.4. Here the system does not 
reorder since at low densities the colloids are far enough 
apart that they interact only weakly and remain in a dis- 
ordered state. The pinning density n p also plays a role in 
determining whether the moving stripe state can form. 
In the inset of Fig. [3Jc) we plot Pq versus n p for a system 
with fixed F p = 1.0 and n = 0.385/ag at F D /F C = 3.0. 
For n p < 0.5/ag, P§ remains low, indicating that the sys- 
tem does not reorder into the stripe state. This result 
shows that a critical amount of quenched disorder, as 
well as a critical strength of disorder, is necessary for the 
moving system to form the stripe state. To show that 
there are no finite size effects for the appearance of the 
four regions, in the inset of Fig.^d) we plot P@ vs system 
size L at Fd/Fc = 3.0 in the moving stripe state. There 
is a slight increase in P§ for the smallest values of L since 
the small systems form fewer stripes, resulting in a lower 
number of topological defects at the interfaces between 
any two stripes. We find that P§ saturates for L > 20, 
well below the size studied throughout this work. 

The moving stripe state forms when the quenched dis- 
order is strong enough to induce plastic deformations 
in the colloid configuration. The effective pinning force 
from the substrate decreases for increasing colloid-colloid 
interactions in a monodisperse system. In a bidisperse 
system, local clustering of colloid species occurs, as seen 
in experiments [3J. Clusters of the weakly charged col- 
loid species B experience a stronger effective pinning 
force than clusters of the more strongly charged colloid 
species A. If the quenched disorder is sufficiently strong 
and dense, the two species will move at different average 
velocities even though they both couple to the external 
drive in the same way. Stripes form as the faster moving 
species begins to collect together. In a co-moving refer- 
ence frame, this state resembles a system of bidisperse 
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FIG. 4: (a) P 6 versus q A /q B for F p = 1.0, n = 0.385/ag, 
n p — 4.0/ao, and Fn/Fc = 3.0. MS: moving smectic region; 
IV: moving stripe (region IV); MC: moving clump region, (b) 
Real space colloid configuration for the moving clump state 
at qA/qB = 11- Small circles: species A; large circles: species 
B. (c) The corresponding structure factor S(k). 



colloids moving in opposite directions, where a laning 
instability has been shown to occur [ll[. In Ref. [ll[, 
the two species were driven in opposite directions with 
no quenched disorder; here, the disparity in velocity be- 
tween the two species is induced entirely by the quenched 
disorder. If the quenched disorder is weak, the colloid 
configuration moves elastically and no velocity difference 
between the two species occurs, so that no stripes can 
form. We have verified that the average velocity distri- 
bution for the two different species (not shown) is differ- 
ent in region III but the same in region II. The stripe 
phases show less ordering than would occur in a system 
with only a single species since topological defects form 
at the boundaries between the phases. 



We now address the effect of the polydipersity on the 
moving regimes by varying qa/qb- A system with the 
same parameters as in Fig. [U[a) and F p = 1.0 acts like 
a monodisperse sample for 1 < qa/<1b < 1-3, where it 
reorders into a moving smectic (MS) state as observed 
in previous simulations [1,0, Hi- F° r 1-4 < Qa/qb < 10 
the system forms a moving stripe state (region IV), and 
for Qa/qb ^ 10 we observe a moving clump (MC) state. 
The onset of these different states appears in Fig. [4ja) 
where we plot Pq versus Qa/qb at Fu/Fc = 3. In the 
moving smectic state, P§ « 0.82, in the moving stripe 
state Pq « 0.7, and in the moving clump state Pq is 
sharply reduced. The real space configuration of the 
moving clump state is illustrated in Fig. QJb), and the 
corresponding structure factor S*(k) in Fig. EJc) has a 
multiple ring structure characteristic of a disordered sys- 
tem containing more than one length scale. Here, the 
different length scales are associated with the triangular 
ordering of the higher charge qA species and the inter- 
clump distance. Similar to the moving stripe state, the 
moving clump state requires sufficient quenched disorder 
to form. As qA/qB increases, the size of the individual 
clumps decreases. 

In summary, we have shown that a bidisperse system 
of repulsively interacting colloids driven over quenched 
disorder can form moving smectic, moving stripe, and 
moving clump states. The moving states are in general 
more ordered than the equilibrium states that form in the 
absence of quenched disorder, in contrast to a monodis- 
perse system. The stripe state arises when the disorder 
is strong enough to induce different average velocities for 
the two species, leading to a laning instability similar 
to that seen for particles driven in opposite directions. 
In our system, the velocity difference is due entirely to 
the quenched disorder, and not to different drives on the 
different colloids. Physical systems in which this model 
could be realized include driven bidisperse colloidal as- 
semblies with charge or magnetic interactions driven over 
optical pinning or a rough wall, electron systems with co- 
existing single and multiple electron bubble states, and 
superconductors with mixtures of Abrikosov and Joseph- 
son vortices. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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